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Abstract. We establish the local uniqueness of steady transonic shock solutions with spherical 
symmetry for the three-dimensional full Euler equations. These transonic shock-fronts are impor- 
tant for understanding transonic shock phenomena in divergent nozzles. From mathematical point 
of view, we show the uniqueness of solutions of a free boundary problem for a multidimensional 
quasilinear system of mixed-composite elliptic-hyperbolic type. To this end, we develop a decom- 
position of the Euler system which works in a general Riemannian manifold, a method to study a 
Venttsel problem of nonclassical nonlocal elliptic operators, and an iteration mapping which pos- 
sesses locally a unique fixed point. The approach reveals an intrinsic structure of the steady Euler 
system and subtle interactions of its elliptic and hyperbolic part. 



1. Introduction 

The study of the Euler equations for compressible fluids is one of the central topics in the 
mathematical fluid dynamics, and the analysis of solutions to the system is of particular interest in 
applications. In particular, in the recent years, important progress has been made in the analysis 
of transonic shock solutions of the steady potential flow equation and the steady Euler system in 
multidimensions (cf . [H El [6l O [HI [9l [151 [211 [22] and the references cited therein) . 

In this paper, we are concerned with the local uniqueness of transonic shock solutions with 
spherical symmetry for the three-dimensional, steady, full Euler system of polytropic gases. Such 
a study not only helps us to understand transonic shock phenomena occurred in divergent nozzles, 
which have many important applications, but also provides new insights for the theory of free 
boundary problems of partial differential equations of composite-mixed elliptic-hyperbolic type. 
This problem can be formulated as a free boundary problem for the Euler system in a spherical 
shell, with the transonic shock-front as a free boundary, which is a graph of a function defined on 
S^ (the unit 2-sphere in M^). Therefore, for such a problem, although there exists a system of 
global Descartes coordinates, it is more convenient to use the local spherical coordinates and the 
terminology of differential geometry (see Appendix A). 
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Let < be two positive constants. A spherical shell in centered at the origin is a 
Riemannian manifold Ai = [r^,r^] x with a metric G = Gijdx'^ dx^ in local spherical 
coordinates (see Appendix A). Its boundary dM is TP U Y} with 

:= {(x°, x) eM:x^ = r\x = (x\ x^) e S^}, i = 0, 1, 

denoting respectively the entry and exit of M. 

Let p, p, and in represent the pressure, density, and entropy of gas flow in the manifold 
respectively. For polytropic gases, p = A{S)p^ with the adiabatic exponent 7 > 1, and the sonic 
speed is 

c = Vipfp- 

Let u be the velocity of fluid flow, which is a vector field in Ai whose integral curves are called 
fluid trajectories in M. Then the steady, full Euler equations for compressible fluids in M are (cf. 
§16.5 in [18]): 

if := div(pn (g) u) + gradp = 0, (1-1) 
ipi := div(pu) = 0, (1.2) 
ip2 ■■= dw{pEu) = 0, (1.3) 
where div and grad are respectively the divergence and gradient operator in A4, and 

1 c2 
E ■= -G{u,u) 



2 ' ' ' 7-1 

We use U = {p,p,u) to denote the state of fluid flow. If U depends only on x^ and u = u^{x^)do, 
then (|l.ip - (|1.3p can be reduced to the following differential equations: 

dnO 2c^u^ 



dxO x0((u0)2 - c2) 
dp 2p{u^f 



dxO x0((n0)2 -c2) 
dp 2pc2(ti0)2 



(1.4) 
(1.5) 

dxO ~ X0((u0)2 -c2)' '■^•^^ 

It has been shown in Yuan [22] that, for equations (|1.4p - (jl.6p . given supersonic data (r^) = 
{Pj^ (r^), Pf^ {r^),u'^ (r^)) on the entry S*^, there exists an interval / such that, if the back pressure 
p^{r^) G I, then there exists a unique S (r*^,r^) so that 

Sb = {ix'^,x\x^) £M : x° =r6,(x\x2) G S^} 

is a transonic shock-front, determined by the Rankine-Hugoniot jump conditions: The flow U^{x^), 
x^ e {r^,ri,), ahead of S}, is supersonic; Uj^{x^),x^ G {rh,r^), behind of Si, is subsonic; and the 
physical entropy condition p'^i^ri,) > pf^irh) holds on S),- We call such a spherical transonic shock 
solution Ub := {U^ ; Sb) to the Euler equations as a background solution. The objective of 
this paper is to study the uniqueness of transonic shock solutions, if it exists, in a neighborhood 
of a class of background solutions under the three-dimensional perturbations of the upcoming 



tin this paper we always use the Einstein summation convention for the Roman indices from to 2 and for the 
Greek indices from 1 to 2. 
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supersonic flow U^^ at the entry S*^. This class of background solutions especially include those 
background solutions satisfying the S-Condition defined in §2. In particular, we prove that, for 
a given back pressure, for sufficiently large upstream pressure and Mach number, the background 
transonic shock-front itself is locally unique. 

The main theorem of this paper is the following: 

Theorem 1.1 (Main Theorem). Let Uh and 7 satisfy the S-condition (see Definition 2.3 in §2j. 
Then there exist Eq and Cq depending only on Ub and 7 such that, if the upcoming supersonic flow 
on S*^ satisfies 

for some a G (0, 1) and there exists a transonic shock solution U = {U^ ,U^; S"^) of (|1.1|) - (!1.3|) in 
M satisfying Property (A) below, then this solution is unique. Here Property (A) consists of the 
three conditions: 

(i) S"^ = e M : x° = ip{x^,x'^),{x^,x'^) € S^} is the shock-front; U~ is the 
supersonic flow ahead of S"^ ; and [/+ is the subsonic flow behind of S"^ . The physical 
entropy condition: p^lgv > P holds on . Moreover, tp G C^'"(S^) satisfies 

IIV' - ^b|lc4."(s2) < C'oe; (1-8) 

(ii) The back pressure is 

p{r\x\x^)=p+{r^)- (1.9) 

(iii) For = {{x^^x^^x"^) e M : ±(x° - V(x\ x^)) > 0, (x\ x^) G S^}, 



U^-Ut\\c^,.^M^)<Coe, (1.10) 



with 

ll^llc"^''^{>!) WiP^ P)\\c'''°'{M) + ll'"llcf'"(>!) ' 
where u is the 1-form corresponding to the vector field u via the Riemannian metric of M, and 
Cr''^ denotes the space of r -forms in Ai (i.e., A^^Ai)) with (7'^'" components in local coordinates 
and the norms are defined in the usual way by partition of unity in the manifold. 



As explained in [10^ [22] , the background solutions coincide with the solutions of the steady quasi- 
one-dimensional model of flows in divergent nozzles. Therefore, the above uniqueness result will 
help to understand transonic shock phenomena in divergent nozzles, as well as the effectiveness of 
the quasi-one-dimensional model [19j . Also see [20j for an explanation of the quasi-one-dimensional 
model from the viewpoint of flows in Riemannian manifolds and [E] for the stability result of 
transonic shock-fronts for the two-dimensional case. 

Apart from the physical implications, the approach by considering the Euler equations in a 
Riemannian manifold is of interest itself in mathematics. We note that some studies have been 
made for conservation laws in general Riemannian manifolds (cf. [H [21 [l7] and the references cited 
therein). This approach via differential geometry reveals some intrinsic structures of the steady 
Euler system, which are valid in general Riemannian manifolds. 

Finally, we remark that the existence and uniqueness of supersonic flow U~ in Ai subject to the 
initial data f/~|so satisfying (jl.7p follow directly from the theory of semi-global classical solutions 
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of the Cauchy problem of quasilinear symmetric hyperbolic systems if £q is sufficiently small (cf. 
[14^ 118]). Furthermore, one can obtain 

where Ci > and Eq > depend solely on {r^) and r^. Thus, we focus on showing the 
uniqueness of and below. Indeed, what we obtain here is much more than this: We design an 
iteration mapping and show that it always has a unique fixed point, and any solution to (jl.ip - (jl.3p 
must be a fixed point of this iteration mapping. Therefore, the solution to (|l.ip - (|1.3p is unique, 
and then Theorem 11.11 is proved. However, we have not known whether the fixed point of the 
iteration mapping is a solution to the original problem (|l.ip ~ (jl.3p . therefore, the existence problem 
for solutions to (ll.ip ^ fll.Sp is still open, though we believe that the ideas and approaches developed 
here will be useful to establish an existence theorem which is out of the scope of this paper. 

For simplicity, we write as U from now on. We emphasize here that the supersonic flow U~ 
is defined in the whole A^, while, by Proposition 9 in [22], can be extended to [rf, — /if,, r^] x 
and still obeys the Euler system, with hi, > depending only on U^{r^). 

The rest of this paper is organized as follows. In §2, we derive some elliptic or transport equations, 
as well as certain equations of exterior differential forms in M and on the shock- front, from the 
Euler equations (jl.ip ~ (jl.3p and the Rankine-Hugoniot jump conditions. Most of the formulas 
obtained here are also valid in general Riemannian manifolds. Based on these decompositions, in 
§3, we present an iteration mapping and establish the existence of a unique fixed point, which 
yields Theorem II. li Some facts and notations of differential geometry are shown and described in 
Appendix A. 

We remark in passing that the analysis and results developed here should be straightforward 
extended to the higher dimensional case, even general Riemannian manifolds for most of them. 

2. Reduction of the Euler system and Rankine-Hugoniot Jump Conditions 

In this section we introduce a reduction of the Euler system and analyze the Rankine-Hugoniot 
jump conditions. 

2.1. The Euler Equations in M. We use d to denote the exterior differential operator and D^uj 
to denote the covariant derivative of a tensor field oo with respect to a vector field u in 7W; while 
is the covariant derivative on S^. CuOJ is the Lie derivative of a; with respect to u in A^. The 
symbol A represents the Laplacian of forms in M., and A' is the Laplacian of forms on S^, which 
are both positive operators (cf. (jA.lOp ). Note that '0 £ C^'"(S^) defines a mapping from to M. 
by (x^, x^) I—)- (^/'(x"^, x^), x^, x^). We use ifj* to denote the pull back of forms and functions induced 
by this mapping; for example, for a function p G ^{M), ip*p = p\si>. The volume 2-form of is 
written as vol, and vol^ is the volume 3- form of Ai. 

In the following, we derive some well-known equations from the Euler system (ll.ip - (11.3p which 
are valid only for flows (cf. [E]). Since these involve differentiations of (ll.ip - ()1.3p . the solution 
of the reduced equations might not be a solution to the Euler system (jl.ip - (jl.3p . One point below 
is to express the relations between these derived equations and the original Euler equations, which 
may be useful in the future to verify that the solution of these derived equations satisfies the Euler 
system indeed. 
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The conservation of mass (jl.2p can be written (equivalently for flows) as 

ipi = DuP + pdiv u = 0. (2.1) 

By the identity 

div(u (8) u) = (div v)u + DyU 
for two vector fields u and v, the momentum equations (jl.ip become 

(/Jo := -(v^ — V^iu) = DuU + -gradp = 0. (2.2) 
P P 



Similarly, the conservation law of energy (jl.3p may be written as 



^3--=-{V2-Eipi) = DuE = 0. (2.3) 
P 

This is exactly the well-known Bernoulli law. 
Since 

^I?„(|up) = G{DuU,u) = G{ipo - -gradp, u) = G{(po,u) - -D^p, 
2 ' ' P P 

we have 



(^3 _ G(^o,n)) =-—(^D4AiS)p'r~') - -D^{A{S)p^) 
p' p'^-K^—l P 

= DuA{S) = G, (2.4) 

i.e., the invariance of entropy along the flow trajectories for flows. 

For a vector field u = u^di in we always use u = u^Gijdx^ to denote its corresponding 1-form 
with respect to the metric G of Ai. Then (12. 2p is equivalent to 



^o = DuU+ — = CuU-di^-^) + ^ = 0, (2.5) 
P 2 p 

since CuU = DuU + d{^-^). This implies 

Cudu = —d(-) A dp + d(fo, (2-6) 
P 

which is a transport equation of vorticity. Moreover, du\gii, the initial value of vorticity on S"^, 
expressed in local spherical coordinates, is 

du\si, = d{^*{v}Gij)) A dx^ - {^*do{u'Gij))dil) A dx^ + {'il^gap^*uP - ^*{^))dx^ A dx" 

+i;^g^pr{^^^)dx'^dx^ (2.7) 

where g = gai3dx°' dx^ is the standard metric of S^. 

Let d* be the codifferential operator in Using the identity d*u = — divu, we obtain 

^-^ = -d*{pu) = 0. (2.8) 

Let (•, •) be the inner product in Ai of forms and let * be the Hodge star operator, which imply 

*vol^ = 1, *1 = vol^, d{*u) = (divn) vol^, a A */? = (a,/3) vol^. 
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Then ([ZH]) is 

(pi = —pd*u + {dp,u) . (2. 
Note that, by the equation of state p = A{S)p'^ , we have 

DuP dp DuA{S) 



{dp, u) = DuP 



2 QJ^ .2 



OA C 



Thus, by setting 



ipi:= — + -^iiP2-Eipi-pG{(po,u)), (2.10) 
p pc^ 

the equation of conservation of mass may be written as 

= _d*n+:^ = divn + :^=0. (2.11) 

IP IP 

By this equation and the definition of the Laplacian of forms A = dd* + d*d, we also have 

= dT {du) + d{^) - d^i. (2.12) 

2.2. The Rankine Hugoniot Jump Conditions on a Shock-Front. A shock-front is a 
hyper-surface in M across which the physical variables have a jump. In our case, it can be expressed 
as the graph of a mapping -.S"^ ^ M.. In local spherical coordinates, we may write 

= {(x°, x\ x^) G 7W : x° = i^{x^,x^), (x\ x^) G S^}. (2.13) 

The normal vector field and the corresponding normal 1-form of S"^ with respect to M. are 

n = {da,4'G'^^dp-G^^do)\s^ and n = {d^^ - dx^)\s^. 

From (ll.ip - (ll.3p . the Rankine-Hugoniot jump conditions (i.e., the R-H conditions) on are 

\G{u,n)pu + pn\\gi, = Q <;=^ \n{u)pu + pn\\g,ii = (2-14) 
\G{u,n)p\\s^=^ ^ \n{u)p\\s.p=0, (2.15) 
\G{u,n)pE\\s^=^ ^ \n{u)pE\\s^ (2.16) 

where [-J denotes the jump of a quantity across S"^ . It is well-known (see lllj) that a piecewise 
G^ state U = {U~ , U^; S^) is a weak entropy solution of (ll.lj) - (ll.3j) if and only if U satisfies these 
equations in Ai^ in the classical sense, the R-H conditions along S'^, as well as the physical entropy 
condition p'^\s^ > P~\ ■ 
Set 

n = no + ni :=n°ao + n"(9„, (2.17) 
n = n° + := u^Gqo dx° + u^'Gafi dx^ . (2.18) 

Then n{u) = dtp{ui) - n°, and (1231) can be separated into 

'4)*\pd'il) - gu^\ = 0, (2.19) 
r\9u''+p\=Q, (2.20) 
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where g = g{U, ip, Dip) := {pu^—dip{pui))\gi, is a function on S^, with dtp £ yl^(S^) being considered 
as a 1-form in Ai; and the expression g{U, ip, Dip) means that g depends on U, ip, and the first-order 
derivatives of ip. Then (j2.15p - (|2.16p become 

r\9\=0, r\E\=0. (2.21) 

Equations ()2.3p and ()2.2ip indicate that E is a constant along the same trajectory even across the 
shock-front. Therefore, we may write E\gi, = Eo{x),x G S^, with Eq{x) a given function depending 
only on the supersonic data at the entry. 

Now, if Ip* \p\ 7^ (which is guaranteed later by the physical entropy condition of the background 
solution), i^Jm yields 

di; = u:= ^^^1^ = florin') + go{U, C/", ^, D^) E ^^(S^) (2.22) 
by using the first equation in (|2.2ip . with 

{puX _7 + l (W)' 



Pb - Pb 



fc2 - (u^^ 



>0, 

St 



and gQ the higher order term defined below (see Definition 12. 2p n. which contains U, U~ ,ip, and Dip 
(the first-order derivatives of ip). We note that equations (|2.20p - (|2.22p are equivalent to (j2.14p - 
()2.16p . Equation ()2.22p also indicates dw = 0. Therefore, we have 

d{iP*u^)=x{U,U-,iP):=-'^. (2.23) 

^0 

Definition 2.1. A constant is called the position of the surface S"^ defined by (j2.13p provided 
that fg2{ip — f^) vol = 0. The function ipP := ip — is called the profile of S^. 

Remark 2.1. The reason why we distinguish the "position" and the "profile" is that they are 
determined by different mechanisms: The "profile" is determined by the R-H conditions, while the 
"position" is determined by the solvability conditions closely related to the conservation of mass. 

Definition 2.2. Let IJ = U'^ — . A higher order term is an expression that contains either 

(i) U" — and its first-order derivatives; 

or 

(ii) the products of ip^, r^ — rb,U, and their derivatives DU, D^U, Dip, D'^ip, and D^ip, where 
D'^u are the /c^'^-derivatives of u in local coordinates. 

Next, we linearize the R-H conditions (|2.20p - (|2.2ip . We write them equivalently as 

Gi {ip* U,rU-) = ^i{rU,iP*U-,DiP), i = 1,2,3, 

with 

Gi = ip*\p{u'>f+P\, ^i = r{\diP{pu''ui)\), (2.24) 

G2 = iP*\pu''\, ^2 = r{\dij{pui)\), (2.25) 

G3 = r\E\, ^3 = 0. (2.26) 



-'^From now on, we always use gi to denote the higher order terms on , and fi to denote the higher order terms 
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As in [15], since Gi{U^ {ri„x),U^ {r},,x)) = for x = {x^^x^) G S^, we have 

d+G,{U+{n,x),U^{n,x)) . {U{i,{x),x) - U+{i;{x),x)) (2.27) 

= { - {d+Gi{U+{i^{x),x),U^{i^{x),x))-d+Gi{U^{rb,x),U^{ri„x))) 

• {U{i;{x),x)-U+{i;{x),x)) 
+d+Gi{U+{i,{x),x),U^{i;{x),x)) . {U{4:{x),x)-U+{i;{x),x)) 

- {Gi{U{tP{x),x), U- {tP{x),x)) - Q(f/+(V(x), x), U- {tP{x),x))) 

- {Gi{U+{i^{x),x), U- Wx),x)) - Gi{U+{ij{x),x), (^(x), x))) + 

-[Gi{U^{ij{x),x),U^{ij{x),x))-Gi{U^{n,x),U,{n,x))] 
=:Ii+IIi, 

where we use "•" as the scalar product of vectors in the phase (Euclidean) space, and d+Gi{U, U^) 
and d-Gi{U, U~) as the gradient of Gi{U, U~) with respect to the variables U and U~ , respectively. 
By the Taylor expansion, the terms in Ij are of higher order. However, 

IIi = ^(p+(rfe)-p-(rb))(V^P + rP-rfe)+0(|^-rb|2), 11,- = 0, j = 2, 3. 
The Landau symbol 0(|^/^p) means the terms of order at least to be two of ^. One can also obtain 
do 



^ 9(Gi,G2,g3) \ 



{U^,U+;St) 



/ 2p+(^0) 



det 



Pi 



1 





((nO 



Thus, ()2.27p is equal to 



Pt 



7-1 

1 





> 0. 



i-X 1 

7-1 p+ 



,.o_.A rp I 



{ptin) -Phin)) 






(ip - rb) + h.o.t.. 



where h.o.t. represents the higher order terms for short. 
We can solve this linear system to obtain 



r {u^) = fii i^pP + rP- rb) + gi {U, , V, D^:), 
r ip) = ^2 (V-^ + rP- rb) + g2{U, , V, D^;), 
r iP) = P3 ir + rP- rb) + g3{U, , i^, D^,) . 



(2.28) 
(2.29) 
(2.30) 
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Using A{S) = pp''^, we also obtain 



riA{S)) = //4 {V + rP- n) + g4{U, U', V, D^), (2.31) 



where 



Ail 



— — — >0, M2 = -7 — — ^( 7 - +c ). (rfe) < 0, 

(7 + 1)^6 (7 + l)r-6(p^(r6))T ^ 



2.3. Restriction of the Conservation of Mass and Momentum on the Shock-Front. We 

now calculate d*{ijj*ui) by restricting the equations of conservation of mass and momentum on the 
shock-front S'^ and obtain a second-order elliptic equation for i/jP. 
In local spherical coordinates, we have 



d*{ru') = -div(V'*KG„^)/^c'^) 

= -V'V*(div« - -^doiVOu^)) - ip*do{ix^fdtp{ui)) 



V G 

= V'r{^ - 'Pi) + rdo{{x'f u') - rOoiix'^fdi^iui)). 
On the other hand, from the conservation of momentum, 

r {dx^ifo)) = r («°5ou°) + r (— ) + r + «°«^ro^) 

P 

-dtp{tp*ui){ij*dou^) -tpg{i;*ui,ilj*ui), 

where we have set ijj*ui := (V'*u°)5ce which is a well-defined vector field on S^, V^*ui (-0*^°) is a 
higher order term, and r% are the Christoffel symbols (see Section A.l). Then we obtain 



d*(^* -til) = 1 + 11 + 111, 
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with 

1 = ^^ (^^0 ^i)=i^^ +h.o.t., 



III = -^{i, + dipirui) (v*5ouO) - v^*uAru'')) 

-rdoiix^diPim)) + ^{V^,^,{rp) - #(rni)(r9op)), 

where, in the term I, we used that ipo = (fQ — {fo)t small so that dx^{(pQ)(-\j — \) is a higher 
order term. 

Note that, for the background solution, II = 0. Then the Taylor expansion and the boundary 
conditions (l2:28D - (f2:3Tl ) yield 

II = /X5 ridop) + ^,r + ir" - n) + 0{\ ir, rP - n, r{U))?). (2.32) 

with 



^5 



<0, 

r=Tb 



'"^ (7 + l)(l-l) ""-'"--^''^''^°- 

where ts = t{r^) = {M^f{rb) G (0,1), and M^{n) = ^""jf^^'' is the Mach number of the flow 
behind the transonic shock-front of the background solution. 

We note that the term III consists of higher order terms. Then we obtain 

+55(C/, U- , ^P, DU, Dij) + i^^r (^^^ - . (2.33) 

Therefore, by (j2.22p . we have 

AV^ + ^uyV'^ = fioiJ-eir^ - rb) + fiofJ-5 i'* dop 

,2 i^fdx'^i'Po) 



with ge = ^055 + d*go and fij = -fiofJ-6 < 0. 

By (j2.33p and the divergence theorem, we choose 



rP-n = --^[ (fi5r{dop) + 95{U,U-,ij,DU,Dij))Yol. (2.35) 
Substituting this into (j2.29p . we obtain 

^p = J_U*p_^^ f ^l;*{dop)vol + griU,U-,i;,Di;)) , (2.36) 
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with fis = — < and 57 = 2 95 vol — §2- Therefore, from (j2.34p . we also obtain an 

equation for the pressure on 5^: 



= gsiU, U',i^, DU, D^, D^U, Z^V, ^'V') + /U2/io^'^* (^^^ - <^i) , (2.37) 
where ^9 = -/io/"2/^5 < and 5-8 = A' 92 + Ai75'2 + + M2C^*5o- 

2.4. An Elliptic Equation for the Pressure in M^. First, we note the following tensor identity: 
div{Duu) - Du{divu) = ClCl{Du (g) Du) + Ric(u,u), 

where Du is the covariant differential of the vector field n in a Riemannian manifold, Ric(-, •) is the 
Ricci curvature tensor, and Cj(r) is the contraction on the upper i and lower j indices of a tensor 
T. In our case, since Ai is flat, Ric(n, u) = 0. From the Euler equations, we have 

div{Duu) - Duidivu) = div cpo - Du(pi + Dui^^^) - div(^^— ^), 

while direct calculation yields that, for = dodo, 



+—Mc^ - {uy){dopf+pdopdo{^'''^' 



\ ^2 

+h{U,DU,D^U), 
97,0 2 1 2 fi/O'l^ 

+f2{U,DU,D''U). 

Now the point is that the above right-hand sides may be expressed as functions of dQp,doP,p, 
ip*A{S), and some higher order terms. 

Indeed, due to the conservation of momentum, 

dou^ = ^dx^iifo) - -^dop + h{U, DU), 

where /s = ■^^^^G{ui,ui) — -^D^-^u^ is a higher order term. From the Bernoulli law, 

doin^ = ^—doic') + ^ + UiU, DU). 

7 — 1 u" 

By the equation of state, 

do{c^) = (7 - l)p~^A{S)ydop + p^A{S)''^doA{S). 
However, the invariance of entropy implies that 

doA{S) = ^ + MU,DU). 
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Therefore, we have 



A{S){x'',x) = rA{S)+ ^{s,x)ds + h{U,i;,DU), (2.38) 



and 

c2 = {rA{S)) ^ + Li((^4) + friU, 4^, DU), 

= 2Eo{x) - :^P^ {rA{S)) ^ + L2(^3, V4) + fsiU, V', DU), 
where Li and L2 do not hivolve the derivatives of (fi, and Li{0) = L2(0, 0) = 0. We then have 

pxl{^div(fo - Du(pi - —Q-'^i +-^3(c^a;°((y9o),'/'3,'/'4)j 

= xl{l- ^^)dlp - A'p + 2x0(2 - ^^)dop 

= ei dip - A'p + 62 dop + 63^ + 64 rp + f9{U, U-,iP, DU, Dip, D^U) (2.39) 

by (fM]) and (lOTll . where ^3(0, 0,0) = and = 6i(x°), i = I,-- - ,4, are known functions 
determined by the background solution: 

61 = {xy{i-t)>o, 

2x^ 

62 = ^((l + 27)t2-3t + 4) >0, 

^3 = ^^^^(6-19t-7t2(-2 + 7) + tS(l + 27)+t3(_3 + 27-472)), 
^2 (7 - l)(i - 1)'^ ' 

^''''(^ + (^^-l'*'.((2,-3)tHa-3), 



where t 



(p+(r5))^(l + (7-l)ts) (1-0^ 
t{x^) = (M+)2(x0) G (0,1) and = t(r,) = (^^^?|^)' = (M+)2(rb) G (0,1) with 



-^b ) ~ . Since 61 > 0, (|2.39p is an elliptic equation for p. 



^6 



We also recall that t{x^) is monotonically decreasing for the background solution and satisfies 
the following differential equation (cf. [22]): 

dt 2t 2 + (7 - l)t 



dxO a;0 t - 1 



(2.40) 



2.5. The Normalization of and Reduced Equations. The above equations and boundary 
conditions are obtained in Ai^ for the given subsonic flow U and the shock-front tp satisfying 
Theorem ll.il The computations are relatively easy for the sake of rather simple metric G. However, 
to show the uniqueness of the transonic shock- front and the subsonic flow behind it, we need to set 
up an iteration mapping to find a new ip G K.^: 
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(2.41) 



for a positive constant ao to be specified later, by solving several boundary value problems in Ad^ 
for any ip £ fC„, and then show that there is a fixed point that is unique. Therefore, it is convenient 
to introduce a C^'^-homeomorphism \I' : {x^,x) G i— )• {y^,y) G := [0, 1] x by 



^ — 'ip(x) ' 



1 ™2\ 



(2.42) 



to normalize A^;^ to n. We set 0^' = {j} x S^, j = 0, 1. Then dn = n° U 0^ We wih use i to 
denote the embedding of in il. We also define the metric of Q to be 



which differs from the metric induced by ^' only those terms involving Dip or ip — ri,. Therefore, 
according to (12.18p . we define, for a vector field u in ft, the corresponding = u'^Gapdy^ . In the 
following, di = in fl for short. Then (j2.7p can be written on as 



du\^o = diu\no) + Ei{U,^,DU, Dtp, D"^ 1^)1^0 A di^ 

daP 



+ (r^ - n) ({r^ - n)nga(SU^\no + ) 

+gg{U, ^, DU, Di,) + (r^ - n) V''g» /^"'f°^ 



\dy^ Ady'^ 



dy^ A dy^, 



(2.43) 



where Ei is a 1-form in fl depending smoothly on U, ip, DU, Dip, and D'^ip. Equations (j2.22p - (|2.23p 
on are respectively 



dip = uj:= fioi*{u') + goiU,U',ip,Di{j) G ^^(S^) 
d{fu') = x{U,U~,iP) :=-^. 

1^0 



(2.44) 
(2.45) 



Similarly, (j2.28p - (|2.3ip are transferred to 



i* (nO) =fi,{r + rP-n) + gi{U,U-,i;, Di^) , 
i*{p) = fi2 ir + rP- n) + 92 (f/, U-,iP, DiP), 
i*{p) = fis ir + rP- n) + gsiU, U-,iP, Dip), 
e{A{S)) = fii + rP- n) + UU, U-,iP, Dip), 



(2.46) 
(2.47) 
(2.48) 
(2.49) 
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where Zxi = > 0. In addition, we have 

— rt, 

,2-*l t 1 ^dy°{ipo) 



, V', DU, DtP) + {r' - n)^^^ - (^1 , (2.50) 
^ A'r + f^7r = f^of^6{r''-n) + ^^^e{dop)+ge{U,U',i^,DU,D4^,D^i;) 

+^01^ I ((n -rb) ^ (2.51) 

mB ^ ''''-n = -ir—l (-r^i*{dop) + g5{U,U-,^,DU,D^))Yol, (2.52) 

([236]) ^ ^P = ]-fi*p-^^ [ i*{dop)vol + griU,U~,tl^,Dtl^)), (2.53) 
([237]) ^ A'(rp) + ^7(i*p) + ^^r(aop) =58(f/,f/",V',^f^,^V',^^^^,^V,^V) 

+/i2MoV' -rb)-^--^yj- (2.54) 

(1239]) ^ (ri-V;)Vo(divv'o-^«'^i-7^I f^^^i+L3((ri-rb)dyO(^o),V53,'^4)) 

V (r^ - rb)y^ + n / 

= ei OqP - (r^ - nf/S'p + (r^ - r5)e2 9op + (r^ - rb)^e3 j5 + (r^ - rfe)^e4 i*p 

-h{U,U-,i;,DU,DtP,D^U), (2.55) 

where , differ from fj , gj by some higher order terms due to the facts that G differs from 
(\I/-i)*(7 by the terms involving ov ip — rb and that tl is smah. We also note that = 
ei{{r^ - rb)y^ + r^) in (I2.55p . 



Remark 2.2. An important observation is that, by (12. 4p . if 994 = 0, we may write 

993 = G{ipQ,u) = (r^ - rfe) • n+((r^ - rfe)y° + r^) • dy°{(po) + h.o.t., 

so we may write L3((r^ — rb) dy^{ipQ),ip3,0) in (I2.55P as an expression of dy^{(pQ) with coefficients 
depending only on y^ by adjusting the higher order term /g. 

2.6. The S-Condition. We now state the S-condition assumed in our main theorem. 
Consider the following boundary value problem: 

eiv" + (r^ - rb)e2v' + (r^ - rb)^(e3 - Xn)v = -{r^ - rbfe^, (2.56) 
^(0) = 1, v{l) = 0, (2.57) 

^'(0) = -^!i±i^(ri-r,), (2.58) 

where A„ = n(n + 1), n = 0, 1, 2, • • • , and Cj = ei{t{y^)) with t(y'') = {r^ — rb)y^ + rb are considered 
as functions of y^ on [0, 1]. 

Definition 2.3. A background solution Ub satisfies the S-Condition if, for each n = 0,1,2, ■■■ , 
problem (|2.56p - (|2.58p does not have a solution. 
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The following lemmas show that there exist certain background solutions which satisfy the S- 
Condition. 

Lemma 2.1. For given 7 > 1, let he a background solution determined by the supersonic 
upstream data p'^ {r^), {r^),t^ {r^) = (M^)^(r'^), and the back pressure p'^ {r^) , where M^{r^) is 
the Mach number of the upstream supersonic flow. Then, for given p^{r^) and the back pressure 
p^{r^), whenp^{r^) and (r^) are sufficiently large, Ub satisfies the S-Condition. 

Proof. We divide the proof into three steps. 

Step 1. By an analysis of the background solution in [22], it suffices to show that, if k = r^— > 
and a = t'^{ri,) = (M^)^(rb) are rather small, then Uh satisfies the S-Condition. We will prove by 
contradiction: we will first assume that (j2.56p - (j2.58p has a solution and then lead to a contradiction. 

We note that, once is given, we can solve t{x^) for all G [r'',r"'^] (i.e., it is independent 

of p^(r''), (r'^), and p^if^))- This can be seen from (I2.40p and the Rankine-Hugoniot condition 
of the Mach number: 

+ Tzl] ( ^ + 7^1\ _ (7 + 1) 



2 



Moreover, since p^{r^) is given, we have the estimate: 

1 < < c for any x° G (r,,, r^), 

Pb i'^b) 

with C depending only on p^{r^), r^, and r^. 

Step 2. We choose ctq such that, for t < ctq, 64 (t) < 0. Hence, for any t'^{rb) < ctq (this requires 
t^-(rO) large), by ^M, 

e4(i(y°))<0 foryOG(0,l). 

In addition, if n > 3, we see that 63 — An < and — — ''''^^^"'''^^^ < 0. So, by the Hopf maximum 
principle, we infer a contradiction if v{l) = for a solution v to (|2.56p - (l2.58p . 

For n = 0, 1,2, we will utilize an energy estimate below to obtain a contradiction if k is also 
small. 

Step 3. We first reformulate Let hn = -^^^^^k- For simphcity, we write the 

independent variable as y. Then, by multiplying 

f r e2(t(s)) 
Pn{y) = exp / (2/i„ + K ) ds 

VJo e.i{t{s))' 

to (|2.56p . we see that w = e~''"^f satisfies 

, d(p„(y)f ) 2 , , , , 2. , ^ 

Pn{y) -, — — + K an{y)w{y) = K Pn{y), 

dy 



where 



. X . .2fi^n + P7f K + P7 e.2{t{y)) e3(t(y))-A^ 
an{y)=Pn{y) 2 77r^ + 



pI P9 ei{t{y)) ei{t{y)) 



Ate) = -(P„te))V'-.fiMi, 

ei(t(y)) 
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By a change of the independent variable y ^ Zn'- 

'-y ds 



with z„ 



the above equation becomes 



Pn{s) 



ds 



(2.59) 



where ' is the derivative with respect to a„ = an{y{zn)), and /3„ = j3n{y{zn)). The boundary 
conditions are 



w{Q) = 1, w'{{)) = 0, w{zl) = 0. 
Now, multiplying w to (|2.59p and integrating on [0, z*] yield 



{w') dz = —K I PnWdz + K j UnW dz. 

Jq Jo 

Note here that, since n < 3, and |/3„|, |a„|, and are bounded by a constant C, 

w'{s)ds dz < 2Cyi*W / "{w'^dz, 



(2.60) 



(2.61) 



(3nw{z) dz 
anw{z)'^ dz 



< C 
< C 







w'{s)dsy dz<C^ I '\w'fdz. 



By (I2:60]l . /o''"(w')^ 7^ 0. Then, from ([2:6TD . we have 



(w')^dz < 



1 - C'k? 



On the other hand, 





2 /• 


/ u)' dz 


<4 / 


Jo 


Jo 



„/\2 . 



This reaches a contradiction when k < min{^=,^=}. Note that C" depends only on r^,r^, 
t^(r'^), n, and p^(r-'^). □ 

Remark 2.3. In Lemma |2.H for the case that k = — r;, > is small, we require only f^(ro) to 
be large (see [22]). 



In the following, we provide some other results on the existence of background solutions that sat- 
isfy the S-Condition. For given [r*^, r^], we note that a background solution is determined by the 
five parameters i'y,rb,p^ in), {rh),t{rh)) with t{ri,) = (M+)2(rb), 7 > 1, n £ (r°,ri), p^in) > 

0, pt{n)>o, and tin) e (0,1). 

Lemma 2.2. For given 7 > 1, p^{rb) > 0, and ctq G (0, 1), there exist a set Si C [0, ctq] of at most 
countable infinite points and a set S2 € [r'^,r^] of at most finite points such that the background 
solution determined 6?/ 7 > 1, G {r^,r^) \ S2, P^i^b) > 0, p^{rh) > 0, and t{ri,) G (0,(7o) \ Si 
satisfies the S- Condition. 
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Proof. Let a = t{ri,) £ [0,ao]. Then, by (|2.40p . we know that t = t{x^,a) is analytical with respect 
to a. Also, by the theory of ordinary differential equations, v = v{y^,a) is an analytical function 
of a. Hence, we establish an analytical mapping : [0, (Tq] — M by 

Ua) = v{l,a), n = 0,1,2,- •• (2.62) 

Now, if /^^({O}) has infinite points, then the zeros of /„ have an accumulate point in [0, do], 
which implies /„ = 0, especially, /n(0) = 0. 

Note that t(x°,0) = 0. Thus, in this case, (j2.56p - ()2.58p become 

((r^ - n)y° + nfv" + 8(ri - n){{r^ - n)y^ + n)v + (r^ - nf {12 - \n)v 

= 12(ri-r,)^4r^>0' (2-63) 

i;(0) = 1, (2.64) 

^'(0) = ~ > 0, (2.65) 

v{l) = 0. (2.66) 

For n > 3, (r^ - r6)2(12 - A„) < 0. By the Hopf maximum principle, we infer a contradic- 
tion at = 0. Therefore, in these cases, /,7^({0}) consists of finite points, which implies that 
U^=3/^H{0}) is countable. 

For n = 0, 1,2, we also recognize that the solution v of problem (j2.63p - (|2.65p is an analytical 
function of the parameter ri,. Hence, we have an analytical mapping Qn '■ [r^,r^] — )• M defined by 
gn{Tb) = v{l,rb). We claim that U^^o5n^({0}) only finite points. Indeed, if ^^^({O}) contains 
infinite points, then Qn = 0, especially gn{r^) = 0. However, in this case, (|2.63p ~ (j2.66p are reduced 
to 

^" = 0' (2.67) 



v{0) = 1, v{l) = 0, v'{0) = 0. 
Obviously, there is no solution to this problem. □ 

Lemma 2.3. For given 7 > l,p^(rt) > 0, and t(rfe) G (0,1), there exists a set C [r'',r^] 
of at most countable infinite numbers of points such that the background solution determined by 
7, Tfc G {r^,r^) \ S3, p'^{rb) > 0,p^{rb), and t{rb) satisfies the S-Condition. 

Proof. As in the above proof, consider v{l) as an analytical function of for each n. If the pre- 
image of the zero has infinite points for some n, then the function is identically zero and, by (j2.67p . 
there is a contradiction. □ 



Lemma 2.4. There exists G (0,1) such that, for given 7 > l^p^iri,) > 0,Pb{fb) > 0, and 
K''^h) = {M^y{ri)) G (0,cTo); there exists r* G (r'^,r-'^) so that the background solution, determined 
6?/ 7 > l,r5 G {r^:,r^),p'^ (riy), p'^ (rb), and t{rb), satisfies the S-Condition. 

Proof. We choose ctq such that, for t < ctq, ei{t) < 0. Hence, for any t{rb) < ctq, 

e4(t(y°))<0 foryOG(0,l) 
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by (j2.40p . In addition, if n > 3, we see that 

\ /A (r^ - rfc)(An + /i7) 
63 - A„ < 0, < 0. 

Thus, by the Hopf maximum principle, we infer a contradiction if v{l) = 0. 

Now, for n = 0,1,2, consider K = U^^Qg^^{{0}) which have at most finite points. Note that 
^ K. Thus, let r^: = sup K < (if K = let r^, = r^). Then the lemma is proved. □ 

This lemma improves somewhat the results of Lemma 12. H but we do not have an estimate of 
— r=K here as that in Lemma |2.1[ 

3. An Iteration Mapping and Decomposition of the Euler System 

In this section we set up an iteration mapping and show that it has a unique fixed point. By 
the derivations in §2, it is easy to see that any transonic shock solution to (|l.ip - (|1.3p satisfying 
those requirements in Theorem 11.11 must be a fixed point of the iteration mapping, which implies 
its uniqueness claimed in Theorem ll.il Another motivation to introduce the iteration mapping is 
for constructing approximate solutions to show the existence of global transonic shock solutions, 
which require further exploration. 

3.1. The Iteration Set. For given tp S /Co-j its position and profile -0^ satisfy 

||'0^||(74,«(s2) < 2cr, \r^ - rf,| < a. 

We solve the candidate subsonic fiow in := {(j;°,j;) e M : x'^ > 'il){x),x G S^}. By ()2.42p . we 
write the set of possible variations of the subsonic flows as 

Os := |?7 = (p,p,m) : ||(p,p)||c.3,«(n) + ll^llc;f"(n) ^ ^ '^o} (3.1) 

with constant 6o to be chosen later. Given U~ satisfying (jl.lip . for any ^ £ fC„ and U G Os, we 
construct a mapping /Cq- x O5 — )• IC^ x Os denoted as 

by the following iteration process. Then we show that T has a unique fixed point in /Co- x Os- 

3.2. A Nonlocal Venttsel Problem for the Candidate Pressure in Q,. We first choose Eq, 
(To, and 60 small enough such that the formulations in §2 valid. For any ijj G /Co-, U £ Os, and U~ 
satisfying (jl.lip . we may express the higher order terms fi and gi in terms of U = + U. 

By (jl.9p and (]2.54p -( r2.55p . we solve p from the following linear nonlocal Venttsel problem: 

eid^p - (r^ - nf/S!p + (r^ - n)e2dQp + (r^ - nfesp + (r^ - nfe^ p\qo 

= h{U,U- ,i;,DU,Di;,D'^U) in il, (3.2) 

p = ou Q}, (3.3) 

/S!{p\^n) + -^dop\n^ + ^l^py^ = UU,U- ,i^,DU,Di^,D^U,D^i,,D^i,) on 0°. (3.4) 

Thanks to Theorem 1.5 in [16J for the Venttsel problem (note that p^^^ < 0) and Theorem 6.6 in 
|13j for the Dirichlet problem, with the aid of a standard higher regularity argument as in Theorem 
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6.19 of |13) . by considering (r^ — 7'f,)^e4 j5|f^o as a nonhomogeneous term and using interpolation 
inequalities, for the solution p G C^'"(0), we have the apriori Schauder estimate 

||p|lc3,a(Q) < C'2( ||p|lco(n) + ||/9||ci,<.(Q) + Il58|lci."(n0)) (3-5) 

with constant C2 depending only on Ub. 

Next, by Lemma [A.4i let Un,m{y) be the eigenfunctions of A' on with respect to the eigenvalues 
A„ = n(n + 1) > 0, n = 0, 1, 2, • • • . Then 

00 2n+l 

P = X] XI '^n,m{y^)Un,m{y)- 
n=0 m=l 

For /g = and = 0, each Vn,m satisfies the nonlocal differential equation: 

ei^n.m + ('■^ - 'rb)(i2v'n^m + (^^ " 'ri,f{ez - \n)Vn,m + {r^ " nf (iiVn,m.{Q) = 0, (3.6) 
Vn,m{'^) = 0, (3.7) 
^J^■9v'n,m{^) + {r^ - n){K + l^l)Vn,m{^) = 0. (3.8) 

First, if Vn,m(0) = 0, (j3.8p says that ^(0) = 0. Hence, by uniqueness of solutions of the Cauchy 
problem of differential equations, we infer Vn,m = 0. 

Now, if Un,m(0) 7^ 0, then, by considering as the unknown, we see that it solves 

^l^n,m + (^^ - n)(i2v'n^^ + {r^ - nfies - Xn)Vn,m = -(^"^ " nfe^, (3.9) 
Vn,m{0) = 1, Vn,M = 0, (3.10) 

V'n,mi0) = -^^^{r'-n). (3.11) 
/^9 

The S'-condition in §2.6 guarantees the nonexistence of a solution to this problem. Thus, Wn,m(0) = 
0, which implies p = 0. 

Therefore, the S-Condition implies the uniqueness of solutions of the Venttsel problem (j3.2p - 
(j3.4p . Then, by p.Sp and a standard argument of contradiction based on the compactness (cf. 
Lemma 9.17 in |13j). we have the apriori estimate: 

||P||C3,"(Q) < C'2(^ ||/9||c>i,«(f^) + ||ff8|lci.«{Q0)) (3.12) 

for any C'^'" solution of the above nonlocal Venttsel problem. Then, by the method of continuity as 
carried out in [16], we see that this problem has a unique solution p G C^'"(ri) satisfying estimate 

3.3. Update of the Candidate Free Boundary. Once we get p, according to (|2.52p - (j2.53p . we 
may obtain the new profile of the free boundary ip^ and the position by 

V^^ = — (plno-^^ / dopyYol + g7{U,U-,i;,Di;)), (3.13) 
fP-n = --^ f (^^dop\no+95{U,U~,i;,DU,Di;))Yol. (3.14) 
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In fact, one may show that Jg2 V'^vol = by (j3.13p and (j3.4p : However, we do not need this in 
the process. We need to improve the regularity of ij:. By (|3.4p and (|3.13p . satisfies this elhptic 
equation on : 

A'r + fi^r = /ioM6(r^ - n) + 5op|no + UU, U-,i;, DU, Dtp, D^i;). (3.15) 

The right-hand side belongs to C^'"(S^). Thus, by Theorem 6.19 in P^, ip = ifjP -\- obviously 
obeys the estimate: 

IIV' - ?'f.|lc4>"(S2) < C'2( ||(55,57)|lc0(C0) + ||/9|lci."(f7) + H^S + 11^6 Ilc2>-(n0) ) (3-16) 

with the aid of ([3T2|) . 

3.4. Solving the Candidate Velocity on ^l^. Now we need to solve u on 0^. To this end, by 
(|2.45p and (|2.50p . we reformulate this problem as 

d{u \qo) = x{U,U ,V):= G C2' (S ), (3.17) 

d*{u'\no) = -j^dop\no+l^S-rb)+MU,U-,i^,DU,D4;) E ^^'"(S^), (3.18) 

where d and d* are respectively the exterior differential and codifferential operator of forms on S^. 
By integrating (j3.4p on S^, we see that the integration of the right-hand side of (j3.18p on also 
vanishes. Then, by Lemma lA.H there exists a unique solution n^l^o € C^'^i^S"^) with the estimate: 

P^bojl^^a.a^ga) < C2(^ ll5o|lc;3."(QO) + || (^5, 56, 5?) Ilc2,<»(n0) + ||/9||ci>'»(Q) + ll^sllcL^CQO) ) • (3.19) 

Therefore, combining this with (|2.46p . the velocity on the candidate free boundary is obtained. 

3.5. Solving the Candidate Entropy, Density, and Vorticity in 0,. We note that the entropy 
can be solved according to ()2.4p and (j2.49p by the following Cauchy problem of the linear transport 
equations: 

DuMS) = in n, (3.20) 

A{S) = f^4{lp-rb)+gi{U,U-,^p,D^p) onQ^. (3.21) 

Note that A{S)^{x^) = A{S)'^{rb) is a constant. By the theory of ordinary differential equations, 
since u is close to {u'^)'^do in C^''^(Q), the trajectories of u still fill 0. We also have the estimate: 

P(5')|b3,a(f^) < C2(^\\Tp - rb||c3,a(S2) + 1 1 54 1 1 (7^," (f^O ) ) " 

Hence, we may solve the candidate density p + by the state function p = A{S)p'^ . Then 



|p||(73.«(Q) <C2{ \\p\\c3.<y(^n) + 



A{S) 



For du = du = du — du^ , by subtracting the background solution from (j2.6p . we formulate a 
linear transport equation: 



1 ^ , „ dp A dp, 
.-)Adp+^^^ 

Pb (P, 



Cudu = -di—) Adp+ , ,7 + fwiU, DU). (3.22) 
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According to (I2.43p . the initial value on Q,^ can be taken as 

M\no = dmno) + Ei{U,7p,DU,D7p,D^iP)\^o Adip 

+ (,1 _ r,)({r' - n)ng^pu%o - )dy' A dy" 

+UU,i^,DU,Di^). (3.23) 
Note that u\qo has been obtained in §3.4: 

u\no = {i\o + (ri - rb)/ii(V' + - n)dy^ + gw{U, U~ Di,). (3.24) 
Therefore, we may solve du from ()3.22p - ()3.23p and obtain the estimate (see Lemma lA.SP : 

\\'^^\\cl'°'(n) ^ C'2( ||/io||c2."(n) + Il59llc2">-(f70) + Il(/5,p)|lc3,-(c) 

+ 11^^ - ''fe||c3."{S2) + ||-U%o||^3,a^g2) + 1 1 5l0 1 1 C3." (^0) ) • 

3.6. Solving the Lower Regular Candidate Velocity in and on il^. From (j2.5p . by sub- 
tracting the background solution, we formulate a transport equation of the velocity u va. (to 
distinguish the lower regular velocity obtained here from the candidate velocity in the next subsec- 
tion, we write u as ui): 



dp_ pdph 
Pt iPt? 



l^u+^i = -^ + T^ + hiiU,^,DU,D^). (3.25) 



Here we used that J~-uiU^ — d{u^ ,ui) = 0. With the Cauchy data ui\^o as the right-hand side of 
(j3.24p . we may uniquely solve ui, particularly its restriction on 0^, i.e., Mil^i. We also have an 
estimate: 

- ||/ll||c-2,a(^) + ||(p,/5)||c3,a(f^) + ||^^^lno||c3.«(S2) 

+ 11^ - ^6llc3.«(s2) + ll5io|lc3'«(nO) )• (3-26) 

3.7. Solving the Candidate Velocity in VL. Note that u; obtained in the above step is only in 
C^'"; it is not our desired candidate velocity. In fact, we will solve the velocity u by the following 
elliptic equation motivated by (I2.12p : 

= d{Uu^^) + d\du) + d{^-^^p)+h2{U,U-,^l:,DU,D^U). (3.27) 

We impose the Dirichlet condition (|3.24p on il^ and the Neumann condition on according to 
(13:25]) by 

D^+U = D^+tii - C^+iii + ( - ^ + ■^P + fii{U,i^,DU,Di^))\^,. (3.28) 



Note that, in local spherical coordinates, 

1 

X' 
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SO it does not contain the derivatives of ■u/. Therefore, we may solve u (i.e., u) in Q by Lemma lA. 21 
to obtain 

ll^llc^"(n) ^ ^^{\\M\c^,-'^(n) + \\M\c^,'"{n) + \\(P^P)\\c^.^in) + I|/i2|lcj-(n) + ||/ii|lc2.«(ni) ) • 



3.8. Well-Definedness of the Mapping T : IC^ x O5 — )■ /Ccr x O^. First we notice that, for any 
ip G JCa, U G C/^ + O^, and satisfying (ll.llj) . it is clear from the definition that a higher order 
term / satisfies 

11/11 <C2(e + 52 + 0. (3.29) 
Then, combining this with the estimates in §3.2-§3.7 yields 

I|f^lb3,-(n) + IIV' - n\\c^,.(^s2) < 02(6^ + + e). (3.30) 
Now we choose Cq = 2C2 and £0 < -^hj- Then, for 6 = a = CqE, the above estimate shows that 
IJ £ Os and £ K.^- 

3.9. Contraction of the Mapping. Now, for z = 1, 2, choose arbitrarily V'^*) G /Co- and f/'^*) G O5, 
and let f/^*^ and ■f/'^*^ be obtained by the above process correspondingly. Then 

||[/« - < C,e{\\U^^~^ - ?7(2)||^,,.(^) + - V(')|b3,.(s2)), (3.31) 

IIV^W - V'^')||c3,.(s2) < C3e(||?7« - f/(2)||c2,.(n) + 11^^'^ - ^(')||c3,<^(s2)). (3.32) 

This can be achieved by employing the equations of ^7*^^^ — U^'^'> and the estimates of higher order 
terms as sketched below. 

First, for any higher order term /, when 5 = a = CqE, ip^'^^ G /Co-, and f/*-*^ G Os, we have 

\\f{^('\u('\u"^'\...)-f{i;('\u('\u~^'\...)\U 

where f/^'^ := U^'^'^ + Ufj' , U~^^^ = ip^"^^* {U~), and || • is the corresponding norm for / when 
ip G C^'" and U G C^'". As two examples, we have 

(a) . By the mean value theorem and (jl.lip . 

\\{^pWy{U- - [/,-) - (^(2))*(f;- _ 

< C\\D{U- - ?7,-)||c2,.(A,)||V(i) - V'^2)||c2,.(s2) < CellV^') - V(2)|b3,.(s2); (3.33) 

(b) . For f^'^ = h{U^\U-'^'\i}^\DU^\D'^U'^^), 

< Ce(||V« - V'(2)||^3,,(g,) + _ (3.34) 

by the definition of higher order terms. 

Next, consider the equations of C/^^^ — ^7*^^^ The right-hand sides of the elliptic equations of 
_ p{2)^ such as f^^ — /g^^ (cf. ()3.2p ) and ^g^'' — ^g^^ (cf. (j3.4p ). are in C". Therefore, we can 



f (1) _ f (2) 
J9 J9 
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obtain a C^'"-estimate of U^^'^ — U^'^\ rather than a C^'"-estimate. The reason is that the loss of 
derivative occurs in solving the transport equations. For example, from ()3.20p and ()3.2ip . we have 

D^,,,{A{si^^ - Msf^) = -D^,,,_^,,,A{sf^ in (3.35) 

(1) -(2) 

A{S) -A{S) 

= ^4 (V^^^) - V'^'^ ) + (54 , U- , , Di;^'^ ) - 54(C/(') , U- , V^') , D^l;^^'^)) on 0°. (3.36) 

Note that the right-hand side of (|3.35p is only in C^'"(0). 

We omit the details of deriving estimates p.3ip - (j3.32p . since the process is similar to those in 
§3.2-§3.7 except the two points explained above. 

Then the mapping T has a unique fixed point if Sq is small enough by a simple generalized 
Banach fixed point theorem. In particular, the uniqueness of the fixed point implies the transonic 
shock solution of ()l.ip - ()1.3p satisfying the requirements in Theorem 1 1.1 1 is unique, as claimed there. 
We also note that, although T always has a fixed point as we proved, it is not clear yet whether this 
fixed point is a solution to ()l.ip - ()1.3p : therefore, in order to obtain the existence result as assumed 
in Theorem II. H it requires further work, which is out of scope of this paper. 



Appendix A. Some Notations and Facts of Differential Geometry 

In this appendix, we present some notations in differential geometry and some basic facts used 
above for self-containedness. 

A.l. The Metric of Ai in Local Spherical Coordinates. In local spherical coordinates, for 
<C < ,0 < < TT, — vr < < vr, the standard Euclidean metric of Ai can be written as 

G = Gijdx' (g) dx^ = dx^ dx^ + {x^ fdx^ (g) dx^ + (x° sinx^)^dx^ (g) dx'^. 

Hence, ^/G := sj det(Gij) = (x")^ sin X . For the Christoffel symbols, since T*-^ — '^fcj' '-'■'^^y the 
following are nonzero: 

T-O _ _™0 pO _ _„0/ ■ 1x2 pi _-pl _ J_ 
^ 11 — ^ ) ^ 22 — ^ Ism X ) , i 01 — i 10 — , 

1^22 = ~ sin x}' cos x}' , ro2 = r|o = ^ , = ^21 = cot x^ . 

We also use (G*-') to denote the inverse of the matrix (Gij), and |np = G{u,u) = GijU^u^ . 
In local spherical coordinates, we write the standard metric of S'^ as 

9 = Saisdx'^ (g dx^ := dx^ (g dx^ + {sin x^)"^ dx"^ (g dx"^. 

Therefore, we have 

/ Vg 

^ := ^det{gij) = sinx^ = 
The nonzero Christoffel symbols are 

722 = — sinx""^ cosx^, 712 = 721 = cot x^. (A.l) 
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A. 2. Some Lemmas. The following results are used in the text. 
Lemma A.l. There exists a unique u G A^{S'^) that solves 

duj = X, d*uj = ip, (A.2) 
i/x G ^^(S^) and ip G ^°(S^) satisfy /s2 X = and /ga ^vol = 0. 

Proof. Since the first Betti number 6i of is (i.e., 6i = 0), by the Hodge theorem, we can 
uniquely solve to via 

Au = d*x + di^ e A\S^). (A.3) 

It suffices to show that ()A.2p holds. 
First, since Ad = dA, we have 

A{duj -x) = dAuj - {dd*x + d*dx) = 0. (A.4) 

In addition, Jga (doj — x) = by the Stokes theorem and the assumption. Since the second Betti 
number 62 of is 1, by the Hodge theorem, the space of harmonic 2-forms on is one- 
dimensional. Note that vol G T-l^ because d* = — * d* and *vol = 1. Therefore, dcj — x = 0. 
Similarly, we have 

A{d*uj-Tp) = d*Auj- {dd*7p + d*dilj) =0 (A.5) 

due to Ad* =d*A, and 

{d*uj - ip) vol = 0, 



by the divergence theorem and the assumption. Note that the zero-th Betti number 60 of is 
1; according to the Hodge theorem, the space T-L^ of harmonic functions on is one-dimensional. 
One easily sees that 1 G Ti.^ . Therefore, d*u: — -0 = as desired. □ 

Lemma A.2. Let the two disjoint components of the boundary dA4 be Ai^ and A4^ , k = 2,3, and 
a G (0,1). Assume that (3 G Cf~^'"(A^), ujq is a C^'°^ 1-form, and uji is a C^^^'°' 1-form in A4. 
Then there exists a unique C^'" 1-form co that solves the following problem: 

Au = 13 in M, (A.6) 
ulj^o = ujqIj^o, (A. 7) 

Moreover, 

Proof. For C M'^, we choose the spherical coordinates: 

= (ri - r,)yO + r,, = y'^, a = 1,2. 
Then we use the standard global Descartes coordinates: 



U ,z ,z 
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Let Lo = Uidz''. By the Weizenbock formula, 

2 

Aw = - ^{dl,LOj)dz^ (A.IO) 

i=0 

holds globahy (cf. [IS])- Therefore, ()A.6p - (]A.8P represent three decoupled boundary value problems 
of the Poisson equations. The uniqueness, existence, and estimates of the solution are then clear. □ 

The following result follows from Lemma 4.6 in |5j. It particularly implies that the norm of a 
smooth function in a manifold 1^2 is equivalent to the norm of its pull back in another manifold 
which is homeomorphic to 

Lemma A. 3. Let Q.i and Q.2 be two open sets in R" and u G C^'"(r22). Let k be a positive integer 
and a G (0, 1). Let $ : ili satisfy ^ £ C'''°'{ni;0,2)- Then uo^ £ C'''°(f]i) and satisfies 

\\u o ^Ilcfc."(ni) < Chllcfe."(n2)' (ATI) 

where C = C{n, ||$||cfc."{ni;n2))- 

Lemma A. 4 ([3j). The eigenvalues of the Hodge Laplacian on are A„ = n{n+l),n = 0,1,2, ■■■ , 
and there are 2n + 1 linear independent eigenf unctions {tin,m}m=i,--- ,2n+i corresponding to A„, so 
that the eigenfunctions {'Un,m}neNu{o},m=i,--- ,2n+i o,"!^^ smooth and form a complete unit orthogonal 
basis o/L2 (s2). 

A. 3. On the Transport Equations Involving Lie Derivatives in Manifolds. In Section 3.5, 
we need solve the differential forms from the Cauchy problems of the transport equations involving 
Lie derivatives. Here we present the basic theorem with a proof. 

Let M be an n-dimensional closed C°° differentiable manifold, Ai = [0, T] x M, X a C'^^^-vector 
field in Ai which is transverse to L* = {t} x M for t S [0,T], and / a C*''-function in (A; is a 
nonnegative integer). Without loss of generality, we assume that X points to the interior of A4 
when restricted on F^. We wish to solve a r-form uj (r > 1) in A4 which satisfies the following 
problem: 

Cxu} + fu = e in M, (A.12) 
w = cjO on r°. (A.13) 

Here C is the Lie derivative in Ai, d is & given r-form in A4, and is a given point- wise defined 
r-form of class on T^. 

We have the following existence and uniqueness results: 

Lemma A. 5. Under the above assumptions, there is a unique r-form u in M that solves ()A.12p - 
(jA.lSp . In addition, there holds 

W^WciM) - ^{ PWc'iM) + ll'^°llcfc(rO) )' (A-14) 
with a positive constant C depending only on \\f\\ck(^j^-^ o^nd \\X\\^k+i(^j^-j ■ 

For the proof, we first get familiar what (|A.12p stands for in a local coordinate chart. 
Let -E be a local coordinate chart of M. Then E = [0,T] x E is a coordinate chart of A^. In 
E, problem (|A.12p - (|A.13P is an initial value problem of the transport equations. To see this, for 
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simplicity, suppose that X = X^d^ + X'^da, with X" = for a = 1, • • • ,n in Ai. Since is a 

global coordinate, X^ = dx^{X) is a C'^-function defined in Ad. By our assumption, X^ is positive 
and bounded away from zero since Ai is compact. 
Suppose that 

UJ = uji^...irdx'^ A ■ ■ ■ A dx'"' = r\ ^ uji^...i^dx'^ A ■ ■ ■ A dx'\ (A.15) 

0<ji <■■■<«,.<« 

e = ei,...i,dx'^ A ■ ■ ■ A dx'^ = rl 6i^...i^dx'^ A ■ ■ ■ A dx'^ . (A.16) 

0<ji<---<v<ri 

Since, for the differential forms a and f3, there hold 

Cxia A (3) = (Cxa) A /3 + a A (£x/3), 
dCxct = Cx{da), and £x/ = Xf (cf. [12j), we obtain 

Cxuj = r\ (x^douJoi^...ir + wo^a-v^oX") A dx^^ A • • • A dx'"- 



0<i2<---<ir<n 

+r\ X^doOJi^...iJx'^ A • • • A dx'^ 

0<il<---<ir<n 

+r\ Y uJoi^-irOaX^dx'' Adx'^ A--- Adx'"- 

(2<---<«r<n 

Y (x°(9owoi2-»r + uioi2-irdoX^^ dx^ A dx'^ A---Adx' 



0<i2<---<ir<n 



0<i2<---<ji<n 



0<ji<---<ir<n ^ ' 



Y (sign(a) a;o.^(,,...i^(^,5,^(,,XO) jdx^i A dx'^ A • • • A dx^^ (A.17) 

Here 'P(r) is the permutation group of {1, • • • ,r}, and sign((T) is the sign of a permutation cr. 
Hence, by dividing X^ from both sides of equation (|A.12p . we have 

If 1 
9oWoi2-i, + (^5oX° + Yo)^Oi2-ir = ^%2-v> <i2 < ■■■ <ir <n; (A. 18) 

and 

/ 



■Jo^h-ir - (^_\)ixo ^ (sign(<^)9i,(i)^%v(2)-v(,.)), < n < • • • < < n. (A.19) 



These are linear transport equations in the x'^-direction. 

We can first solve (|A.18p in E by using the initial data wo«2---j,. Ir" = {'^^)Qi2---ir ^"^^ then substitute 
woia-.-v in the right-hand side of ()A.19P to solve uji^-.-i^ in E with initial data Wijij-'-vIro = {'^^)hi2---ir- 
By antisymmetry of the lower indices, we obtain all the coefficients cj^j.-.j^. Since the quantities in 
()A.12p - ()A.13P are defined globally, oj = a;i^...j,,dx*i A • • • dx*'' we solved is also well defined in Ai. 
Estimate ()A.14p is obvious from these initial value problems of ordinary differential equations. 



LOCAL UNIQUENESS OF TRANSONIC SHOCKS 



27 



For the general case that X ^ X^Oq, equation ()A.12p would be a first-order hyperbolic system 
with ("^^) unknowns. In addition, we can not use the rather simple coordinate charts like E, since 
the characteristic curves (i.e., integral curves of X) may escape E at some t < T by solving the 
initial value problem. 

We note that the Lie derivative behaves well under the homeomorphisms of differentiable mani- 
folds. 

Proposition A.l. Let ^ be a homeomorphism of M. Then 

{'^-^rCxuj = U,x{i'^'^T^)- (A.20) 
This can be shown by using the Cartan formula £x'^ = dix^ + ixdu; and the formulae 

(cf. |12)). where ix(^ is the interior product of oj and X. 

In addition, by Lemma lA.31 the C'^-norm of a differential form in Ai is equivalent under C'^- 
homeomorphism of the manifold. Thus, to prove Lemma IA.5I for the general case, it suffices to 
straighten the vector field X to the form X^do globally by a suitable homeomorphism of Ai. 

Proposition A. 2. For given C''^^ -vector field X in M. which is transverse to {x^^ x M for every 
G [0, r], there is a C'''^^ -homeomorphism ^: A4 ^ A4 such that, for any fixed x^ G [0, T], it is 
also a homeomorphism of {x^} x M, and <1>*X = {^^^)*{dx^{X))do. 

Proof. Step 1. Let E he a coordinate chart of A4 as introduced above and X = X^di with X^ > 0. 
Since X^ is a nonzero function in Ai, we may define a C'^^^-vector field X in 7W by 

X = x'di = do + ^a,. (A.21) 

Step 2. Now X generates a fiow (pt in M by theory of ordinary differential equations since 
Ai is compact: For any P G M, ^{t) = (j)t{P), t £ [0,T], is a curve in At with initial value 
7(0) = (0, P) G ro. We then define <^ : M ^ M hy 

^MP)) = {t,P), P^M. (A.22) 

Note that <t)t{P) G P* since = 1. We now show that <I> is a homeomorphism. 

Step 3. For any Q G it is easy to see that there is uniquely a pair {t,P) with P £ M 
and t G [0, T] such that (ptiP) = Q by solving backward the integral curve of X through Q. So 
<I> is defined for all the points in Ai. In addition, <I> is obviously surjective and injective by the 
uniqueness and existence results of the initial value problem of ordinary differential equations. By 
continuous dependence on the initial data (0, P) and t, we see that $ and are also continuous. 
If X G C^^^, then <I> and <I>~^ are C'^+^-mappings by C'^+^-dependence of solutions of ordinary 
differential equations on t and initial data. This proves that <I> is a C'^'^'^-homeomorphism. 

Step 4- Now, by definition of push-forward mapping (tangent mapping) of vectors, we have 

^>,(X) = ^^X^X) = ($-i)*(XO)$,(X) 

= (c^-i)*(XO)|c^(7(t)) 

= ($-i)*(X°)9o. (A.23) 
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This completes the proof. □ 
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